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Abstract: The optimal state determination (or tomography) is studied for 
a composite system of two qubits when measurements can be performed on 
one of the qubits and interactions of the two qubits can be implemented. 
The goal is to minimize the number of interactions to be implemented. The 
algebraic method used in the paper leads to an extension of the concept of 
mutually unbiased measurements. 

PACS numbers: 03.67.-a, 03.65.Wj, 03.65.Fd 

Key words: State determination, reduced density, unbiased measurement, 
minimal realization, unbiased subalgebra, Pauli matrices. 



1 Introduction 

An n-level quantum system is described by an n-dimensional Hilbert space 7i, or equiva- 
lently by the algebra M n ((D) of the n x n complex matrices. When an orthogonal basis of 
"H is chosen, operators acting on "H correspond to n x n matrices. A positive operator p of 
trace 1 is called state. If we choose and fix an orthonormal basis {ei, e^, ■ ■ ■ , e n }, then a 
state p is determined by the matrix elements pij = (ei\p\ej). Determination of p involves 
n 2 — 1 real parameters, namely, pa (1 < i < n — 1), Re p^ and Im p^ (1 < i < j < n). 

A von Neumann measurement on the system is a family M. = {Pi,P2, . . . ,Pd} of 
pairwise orthogonal projections such that £\ Pi = I. When the measurement M. is 
performed in the state p, the outcome 1 < j < d appear with probability pj = Tr pPi for 
each j EI- Independent measurements on several copies of our quantum system give 
the relative frequencies fj for each outcome j and fj is an estimate of the probability pj. 

Supported by the Hungarian Research Grants OTKA T042710, T063066 and T032662. 



1 



The repeated measurement provides d — 1 degree of freedom concerning the the density 
p, since YljPj = 1- The information we obtained is maximal if d = n which means that 
all the projections Pj are of rank one. p is determined by n 2 — 1 parameters, hence at 
least n + 1 different measurements are to be performed to cover all degrees of freedom. 
Of course, the n+1 different measurements are sufficient in the case when they provide 
"non-overlapping" information. 

The literature of state tomography is very rich, there are several protocols, and the 
efficiency of state reconstruction can be increased if the later measurements depend on 
the outcomes of the former ones [U E] • 

The composite system of two qubits is a 4-level quantum system which is described 
on the space (D 4 = (D 2 ® (D 2 . A state is described by 15 real parameters. Therefore, at 
least 5 kinds of elementary measurements should be made to determine the state of the 
system. 

Denote by A and B the two qubits. Then M 4 (C) = B(H A ) ® B(H B ), where B(H A ) 
and B{Hb) are isomorphic to M 2 ((D). Assume that we can perform measurements only on 
the qubit A. If the total system has the statistical operator p AB , then we can reconstruct 
the reduced density 

(1) rp 

py := TybPab 

after some measurements. In order to get more information, we switch on an interaction 
between the two qubits. If H is the Hamiltonian, then the new state is 

e iH p AB e- iH = W lPAB W* (1) 

after the interaction. (For the sake of simplicity, the interaction is kept for a time unit.) 
The new reduced density is 

p { A ] :=TtbW iPab W*. 

This procedure may be continued by using other interactions and ends with a sequence 
of reduced states p^\p^\ ■ ■ ■ , Pa ■ We want to determine the minimal k such that this 
sequence of reduced densities determines p A s- In other words, we want to minimize the 
number of interactions between the two qubits. It turns out that the minimum number 
is 5. 

Minimal realizations play an important role in systems theory, too [Sj, because they 
represent the state of the system with the minimum possible number of parameters. 
Minimal realizations are known to be jointly controllable and observable for most of the 
known system classes. The above problem of finding the minimum number of reduced 
states can be regarded as a minimal representation problem for a system that consists 
of a pair of coupled qubits. 

2 Algebraic formulation 

Instead of the transformation (0) of the density matrix p A B, we can change the subalgebra 
and we have an equivalent algebraic formulation. The total system is described by the 
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algebra M 4 ((D). We look for subalgebra Ai, A 2 , ■ ■ ■ , Ak such that 

1. Each Aj is algebraically isomorphic to M 2 ((D), 1 < j < k. 

2. The linear span of the subspaces Ai, A 2 , . . . ,Ak is M 4 ((D). 

Given a subalgebra Aj, there is a unitary Wj such that W*AjWj = B(Ha) <S> (Dig. The 
reduced density pj G Aj is the same as the reduction of WjPABW* to the first spin. 
Therefore, instead of the reduction of density after the interaction, we can work with 
the reduced density of pab in Aj. The second condition makes sure that the reduced 
densities in Ai, A2, ■ ■ ■ , Ak determine pab completely. 

The traceless subspace of M4((D) has dimension 15, while the traceless subspace of 
Aj has dimension 3, therefore we need k > 5 to fulfill the requirements. It will turn out 
that k = 5 is possible. 

The algebra M2((D) is linearly spanned by the Pauli matrices: 
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Recall that they satisfy the multiplication rules 

3 

OiOj = SijI + i 2J QjfcO-fc (1 <i,j < 3), (2) 

k=l 

where e^k is the Levi-Civita tensor: 

{0 if 3j, k such that ij = ik, 
1 if the permutation (iii 2 ■ ■ -i n ) is even, 
— 1 if the permutation (iii 2 ■ ■ - i n ) is odd. 

The above rules can be essentially simplified by posing the following two requirements: 

1. aj is a self-adjoint unitary (1 < j < k) and a 3 = —iaia 2 . 

2. a x a 2 + a 2 a x = 0. 

When a triplet (Si, S 2 , S 3 ) satisfies these condition, it will be called a Pauli triplet. For 
such a triplet TrS^ = and Tr SiSj = for i ^ j. The latter relation is interpreted 
as the orthogonality of Si and Sj with respect to the Hilbert-Schmidt inner product 
(A, B) := TtA*B. Furthermore, it can be seen that the two relations above imply (J2J). 

Given a Pauli triplet (5*1, 5*2, S3), the linear mapping defined as 

a i-> I, oi i-> Si, a 2 S 2 , a 3 i-> -iS 1 S 2 

is an algebraic isomorphism between M 2 ((D) and the linear span of the operators /, Si, S 2 
and S3. 

In the algebra M 4 ((D), the elementary tensors cr^ <g) aj form an orthogonal basis (0 < 
hj < 3). All these operators are self-adjoint unitaries and can be chosen to be Si's. 

The next proposition is the main result of the paper. 
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Proposition 1 There are 5 sublalgebras of B(1~La) <8> B(Hb) such that each of them is 
isomorphic to M 2 ((D) and the reduced states determine an arbitrary state pab of the two 
qubits A and B . 

Proof. First we take the following Pauli triplets consisting of elementary tensors: 



{(T 3 <g) <7i, <7i 



5 0-3 


,0-1 


<8> o- 2 } 




n 
u 


1 
1 


n n 
u u 




1 
















1 


1 








1 




o-i, 




3> o"o} 




n 
u 


1 
1 


u 


u 


1 


n 
u 


n 
u 


n 
u 








- 


-1 








-1 





02, 


cr 3 C 


g) 0-2} 




" 





1 " 










1 




1 













1 








{a <g) (J 2 , (7 2 ® (T 3 , <7 2 <g) (7i} 




i 





-i 





i 





-i 








-1 

1 



-1 










1 
10 
10 
10 





1 

1 

-10 





i 








-i 



-1 







-i 

i 
















— i 








i 








— i 








i 


n 
U 


U 


U 








— i 


" 











— i 


i 











n 
U 


i 


U 


U 





— i 





" 


i 




















i 








— i 





" 








— i 








— i 








i 








i 












Together with the identity, each triplet linearly spans a subalgebra Aj (1 < j < 4). 
It is important to observe that all the matrices have vanishing diagonal, moreover the 
matrices are pairwise orthogonal, therefore they are linearly independent. 

If we find another Pauli triplet (Si, S2, S3) such that the diagonals are linearly inde- 
pendent, then we have a fifth algebra A$ such that {Ak : 1 < k < 5} spans linearly 
M 4 ((D). Indeed, if A is any matrix, then we can find T G A$ such that A — T has 
diagonal and this is in the linear hull of {Aj : 1 < j < 4}. It follows that the reduced 
densities in {Aj : 1 < j < 5} determines pab uniquely. 
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Here is an example of the above described triplet: 



1 
2 



1 i 

-1 -1 

-1 1 

-i i 



-1 
i 

— i 



-1 



—l 



—l 



i 1 

1 i 

-i 1 

1 -i 



-1 



These matrices are not elementary tensors (but they are Hadamard matrices JT] up to 
a constant multiple and were found by means of an exhaustive search algorithm on a 
computer) . □ 

Proposition 2 If all the matrices of the Pauli triplet generating the subalgebras Aj are 
of the form io^ <g) cr/ (0 < k, I < 3), then we need at least 6 triplets to span M^C). 

Proof. Assume that a Pauli triplet (Ti,T 2 ,T 3 ) in M±(<\T) is such that every element is 
of the form ±Oj <8> o~j (0 < i, j < 3). 

If Ti = ±<Tj ® Oj and T 2 = ±0^ <S> o~i, then 

=ti2~3 = — '^^{eikmtjlnO'm ® 0~ n ) 

m,n 

jlnO~0 ® O n 
n m 

Since T 3 is self-adjoint but io"j (g> cr^ is not, it follows that exactly one of the relations 
i = k and j = I must hold. At least one of the operators Tj should be of the form ctq <8> o~j 



or cr 



cr . 



We have three operators in the form a^®aj and three in the form Ojtgxxo (1 < J < 3). If 
we have 5 Pauli triplets, then at least one should contain two of the above tensor products 
(up to a sign). One can see that oo <S> Oj and Oi ® o"o cannot be in a triplet, therefore 
a triplet contains two operators in the form oq <g> cr^ or two operators like o~j <S> cr . In 
both cases, the third operator has similar form. Hence one of the operators o"o <8> o~j and 
CTj ® cr appears in two triplets and in this case 5 triplet cannot span the whole space. 

Six subalgebras described in the proposition can be given by the following Pauli 
triplets: 

{<7i (8)0-1,0-!® o 2 , o ® o 3 }, 
{o 2 <g> o 2 , o 2 ® o 3 , o ® oj, 

{O3®03,0 3 (g)0i,0o(g)0 2 }, 

{o 2 (8) o 2 , o 3 ® o 2 , 01 ® o }, 
{o 3 <g) o 3 , Oi ® o 3 , o 2 ® o }, 
{01 (g) Oi, 2 (g) Oi, o 3 <g) o }. 



(3) 



Together with / each triplet linearly spans a subalgebra Aj (1 < j < 6) and the 6 
subalgebras linearly span the whole M 2 (C) ® M 2 ((D). □ 
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3 Generalizations 



Mutually unbiased bases (or measurements) are interesting from many point of view 
jH El EI] and the maximal number of such bases is not completely known ^2] • The above 
discussed setting of state determination is somewhat similar. In this setting we may 
look for essentially orthogonal non-commutativ subalgebras while unbiased elementary 
measurement are given essentially orthogonal maximal Abelian subalgebras, see Prop. 
2.2 of [B]. The next statement is an analogue of Parthasarathy's proposition. 

Proposition 3 Let A\ and A 2 be subalgebras o/M„((D) and assume that they are iso- 
morphic to Mfc((D). Then the following conditions are equivalent: 

(i) If P E Ai and Q G A2 are minimal projections, then TrPQ = n/k 2 . 
(ii) The subspaces A\ (DJ and A2 CJ are orthogonal in M n ((D) 

Proof. It follows from the conditions that n = mk and Tr P = Tr Q = m for the minimal 
projections. Therefore, condition (i) is equivalent to Tr (/ — kP)(I — kQ) = which 
means that (/ — kP) _!_(/ — kQ). Since the subspaces in (ii) are linearly spanned by 
these operators, the proposition follows. □ 

Now we generalize Prop. 2 for n qubits. 

Proposition 4 If all the matrices of the Pauli triplet generating the subalgebras Aj of 
M 2 n ((D) are of the form icr^m ... crfc( n ) (0 < k(i) < 3, 1 < % < n), then we need more 
than (2 2n - l)/3 triplets to span M 2 «(C). 

Proof. First note that a matrix iofcn) ... er fc( - n ) has only real elements or only 
imaginary elements. Among the three matrices of a Pauli triplet (Ti, T 2 , T 3 ), there is one 
imaginary or there are three imaginary matrices. Let N be the number of triplets with 1 
imaginary matrix and M be the number of triplets with 3 imaginary ones. If the N + M 
triplets with identity linearly span the self-adjoint subspace, then 3(iV + M) + 1 > 2 2n . 
Assume that 

N + M = ^—. (4) 

Since the dimension of the subspace of self-adjoint matrices with imaginary elements is 
(2 2 ™ - 2 n )/2, we must have 

N + 3M = . (5) 

One can see that equations (jlj) and (0) do not have integer solution. □ 

We call a family Aii,A4 2 , . . . ,Aid of subalgebras strongly mutually unbiased if the 
conditions in the proposition hold for any pair. The maximal number of strongly mu- 
tually unbiased subalgebras is not know to us even in the simplest case when the large 
algebra is M 4 ((D) and the subalgebras are isomorphic to M 2 ((D). 
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Following [H], we may call a family M. 1) JH2, ■ ■ ■ , M.d of subalgebras weakly mutually 
unbiased if the subspaces M.\ (DJ, M.2 CI, . . . , M. d (DJ are linearly independent. 
We showed that when the large algebra is M 4 ((D) and the subalgebras are isomorphic to 
M 2 ((D), then the maximum number of weakly unbiased subalgebras is 5. 



4 Discussion and conclusions 



The optimal state tomography has been studied for a composite system of two qubits 
when measurements can be performed on one of the qubits and interactions of the two 
qubits can be implemented. Equivalently, we found physically realizable minimal set of 
reduced densities. The transformation described by is realized by a properly designed 
measurement apparatus (see the experimental devices in [9 ). Therefore, an additional, 
but still unsolved problem is to find the unitary to each of the reduced densities that 
transforms the state to be measured to another one belonging to the reduced density 
in question. (Some preliminary results on the effect of the real 4-dimensional rotation 
matrices on the reduced densities is found in [TU].) 

The construction of 5 Pauli triplets of 4 x 4 matrices from the tensor products of Pauli 
matrices contains heuristic steps combined with an exhaustive search for the missing 
final triplet. This makes it practically impossible to generalize the method for higher 
dimensions. Even if the heuristic steps were removed the resulted algorithm would fall 
to the NP-hard category because of the exhaustive search. 

The questions we posed for two qubits can be asked about three (or more qubits). The 
dimension of M 8 ((D) (D/ is 63, so we need at least 21 Pauli triplets to span the whole 
space. In the moment we can construct 22 spanning Pauli triplets by ad hoc method 
(this will be discussed in another publication). 

The problem investigated here motivates the definition of strongly and weakly unbi- 
ased subalgebras. In relation with them, several questions can be raised. 
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